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Abstract 

The nonlinear asymmetric Kelvin-Helmholtz 
stability of the cylindrical interface between the 
vapor and liquid phases of a fluid is studied when 
the phases are enclosed between two cylindri¬ 
cal surfaces coaxial with the interface, and when 
there is mass and heat transfer across the inter¬ 
face. The method of multiple time expansion is 
used for the investigation. The evolution of am¬ 
plitude is shown to be governed by a nonlinear 
first order differential equation. The stability cri¬ 
terion is discussed, and the region of stability is 
displayed graphically. Also investigated in this 
paper is the viscous linear potential flow. 

Keywords Kelvin-Helmholtz stability, Mass 
and heat Transfer, Cylindrical flow. 0 

1. Introduction 

In dealing with flow of two fluids divided by 
an interface, the problem of interfacial stability is 
usually studied with the neglect of heat and mass 
transfer across the interface. However, there are 
situations when the effect of mass and heat trans¬ 
fer across the interface should be taken into ac¬ 
count in stability discussions. For instance, the 


phenomenon of boiling accompanies high heat 
and mass transfer rates which are significant in 
determining the flow field and the stability of the 
system. 

Hsieh [1] presented a simplified formulation 
of interfacial flow problem with mass and heat 
transfer, and studied the problems of Rayleigh- 
Taylor and Kelvin-Helmholtz stability in plane 
geometry. 

The mechanism of heat and mass transfer 
across an interface is important in various indus¬ 
trial applications such as design of many types of 
contacting equipment, e.g., boilers, condensers, 
pipelines, chemical reactors, and nuclear reac¬ 
tors, etc. 

In the nuclear reactor cooling of fuel rods by 
liquid coolants, the geometry of the system in 
many cases is cylindrical. We have, therefore, 
considered the interfacial stability problem of 
a cylindrical flow with mass and heat transfer. 
Nayak and Chakraborty[2] studied the Kelvin- 
Helmholtz stability of the cylindrical interface 
between the vapor and liquid phases of a fluid, 
when there is a mass and heat transfer across 
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the interface, while Elhefnawy[3] studied the ef¬ 
fect of a periodic radial magnetic field on the 
Kelvin-Helmholtz stability of the cylindrical in¬ 
terface between two magnetic fluids when there 
is mass and heat transfer across the interface. 
The analysis of these studies was confined within 
the frame work of linear theory. They both found 
that the dispersion relations are independent of 
the rate of interfacial mass and heat transfer. 
Hsieh[4] found that from the linearized analysis, 
when the vapor region is hotter than the liquid 
region, as is usually so, the effect of mass and 
heat transfer tends to inhibit the growth of the 
instability. Thus for the problem of film boil¬ 
ing, the instability would be reduced yet would 
persist according to linear analysis. 

It is clear that such a uniform model based 
on the linear theory is inadequate to answer the 
question of whether and how the effect of heat 
and mass transfer would stabilize the system, but 
the nonlinear analysis is needed to answer the 
question. 

The purpose of this paper is to investigate the 
Kelvin-Helmholtz asymmetric nonlinear stability 
of cylindrical interface between the vapor and 
liquid phases of a fluid when there is a mass and 
heat transfer across the interface. 


first order theory and the linear dispersion re¬ 
lation are obtained in Sec.3. In Sec .4 we have 
derived second order solutions. In Sec.5 a first 
order nonlinear differential equation is obtained, 
and the situations of the stability and instability 
are summarized. In Sec.6 we investigate linear 
viscous potential flow. In Sec.7 some numerical 
examples are presented. 

2. Formulation of the problem and basic 
equations 0 

We shall use a cylindrical system of coordi¬ 
nates (r, 6 , z) so that in the equilibrium state 
z—axis is the axis of asymmetry of the system. 
The central solid core has a radius a. In the 
equilibrium state the fluid phase ”1”, of den¬ 
sity p^ 1 ), occupies the region a < r < R, and, 
the fluid phase ” 2”, of density , occupies the 
region R < r < b. The inner and outer fluids 
are streaming along the z axis with uniform ve¬ 
locities U\ and U 2 , respectively. The temper¬ 
atures at r = a, r = R, and r = b are taken 
as Ti,To, and T 2 respectively. The bounding 
surfaces r = a, and r = b are taken as rigid. 
The interface, after a disturbance, is given by 
the equation 

F(r,z,t) = r — R — p(9,z,t) = 0, (2.1) 


The nonlinear problem of Rayleigh-Taylor in¬ 
stability of a system in a cylindrical geometry is, 
however, studied by the present author in (Lee [5- 
6 ]). 

The multiple time scale method is used to ob¬ 
tain a first order nonlinear differential equation, 
from which conditions for the stability and in¬ 
stability are determined. 

In more recent years,Awashi, Asthana and 
Zuddin[7] considered a problem in which a vis¬ 
cous potential flow theory is used to study the 
nonlinear Kelvin-Helmholtz instability of the in¬ 
terface between two viscous ,incompressible and 
thermally conducting fluids. 

The basic equations with the accompanying 
boundary conditions are given in Sec.2. The 


where p is the perturbation in radius of the inter¬ 
face from its equilibrium value R, and for which 
the outward normal vector is written as 


VF r (l dp 

11 “ W¥\ - \ 1+ VrM 


+ {*L 

1 dz 


2n -1/2 


(2 - 2) 

we assume that fluid velocity is irrotational in 
the region so that velocity potentials are and 
for fluid phases 1 and 2. In each fluid phase 

V 2 ^)=0. (j = 1,2) (2.3) 

The solutions for = 1,2) have to satisfy 

the boundary conditions. The relevant bound¬ 
ary conditions for our configuration are 

(i) On the rigid boundaries r = a and r = b: 
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The normal field velocities vanish on both cen¬ 
tral solid core and the outer bounding surface. 

d<jP^ 

dr 

d(j) 

dr 

(ii) On the interface r = R + p(9, z, t ): 

(1) The conservation of mass across the inter¬ 
face: 


= 0 

on r = a, 

(2.4) 

= 0 

on r = b, 

(2.5) 


V 

f 90 

\ dr 


7 Ml 

r + V0- VFj 

= 0 , 

dp 

1 dp dcf) 

dp d(f> 

dt 

r dO d6 

dz dz 


OT P\ “oT Ei7 — T7a/i~o77 ~ 

( 2 . 6 ) 

where [ hj represents the difference in a quantity 
as we cross the interface, i.e., |[ hj = h^ 
where superscripts refer to upper and lower flu¬ 
ids, respectively. 

(2) The interfacial condition for energy is 


Lp(1) + W(1) ' VF ) = 5(??) ’ (2 ‘ 7) 
where L is the latent heat released when the fluid 
is transformed from phase 1 to phase 2. Phys¬ 
ically, the left-hand side of (2.7) represents the 
latent heat released during the phase transforma¬ 
tion, while S(rf) on the right-hand side of (2.7) 
represents the net heat flux, so that the energy 
will be conserved. 


In the equilibrium state, the heat fluxes in 
the direction of r increasing in the fluid phase 
1 and 2 are —K\(T\ — Tq) / R\og{a / R) and 
—K 2 (To — T 2 )/Rlog(R/b), where K\ and K 2 are 
the heat conductivities of the two fluids. As in 
Hsieh(1978), we denote 

SM = K ^< r ° ~ Fl 

(R + 77 )(log b - log (R + 77 )) 


Ki(T\ - To) 


(R + 77 ) (log(i? + 77 ) - log a) ’ 
and we expand it about r = R by Taylor’s ex¬ 
pansion, such as 


S(p) = S(0) + pS'(0) + ±p 2 S"(0) 
and we take 5(0) = 0, so that 


(2.9) 


K 2 (T 0 -T 2 ) 

R\og(b/R) 


Ki{T\ — T 0 ) 
Rlog(R/a) 


G{ say), ( 2 . 10 ) 


indicating that in equilibrium state the heat 
fluxes are equal across the interface in the two 
fluids. 


From (2.1), (2.7), and (2.9), we have 
(jj ( dcj)9) drj 1 drj d^ 1 ) dp dcj)9) \ 

\ dr dt r d6 dO dz dz ) 

= a(p + a 2 p 2 + «3?? 3 ), (2.11) 


P 


where 


a = 


G\og(b/a) 


LR\og{b / R) log (R/a) 


a 2 = 77 - a + 


R 


1 


1 


OL 3 = 


R 2 


2 

11 

~6 


log (b/R) log (R/a) 

2 log(i? 2 /a&) 
log (b/R) log {R/a) 


+ 


log 3 (b/R) + log 3 (R/a) 
{log(V R) log(i?/a )} 2 log( 6 /a)_ ' 


(3) The conservation of momentum balance, 
by taking into account the mass transfer across 
the interface, is 

p«(V</> (1) • VF) + V</> (1) ■ Vf) 

= • VF) + V^ 2) ■ VF^j 

+(P 2 ~ Pi + tV • n)|VF| 2 , (2.12) 

where p is the pressure and a is the surface ten¬ 
sion coefficient, respectively . 0 By eliminating the 
pressure by Bernoulli’s equation we can rewrite 
the above condition ( 2 . 12 ) as 
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f deb 1 

'’{at + 2 


d(j) 

dr 


+ 


1 d(j) 

r 89 


+ 


8(j) 

dz 


- 1 + 


1 drj 
r 86 


+ 


drj 

dz 


( dej) dr/ 1 d<t> drj dcj)\ (dr] d(j> dr] 1 dc/r dr] d(j> 

\ dz dz r 2 86 86 dr J 1 dt dz dz r 2 86 86 dr 


a 


(R + rj)\X7F\ 


1 + 1-52 

1 r 86 


|VFP 


<7 


|VF? 


d 2 V ( , (1 dv 

dz 2 \ \r 89 


2 drj 8 2 r] dr] 1 d 2 r] 
r 2 86 dddz dz r 2 dd 2 


1 + I §2 

dz 


(2.13) 


When the interface is perturbed from the equi¬ 
librium r] = 0 to r] = A exp [i(kz + rnO — cut)], the 
dispersion relation for the linearized problem is 


D(uj, k , m) = aocu 2 + (ai + ibi)uj + <12 + ib 2 = 0, 

(2.14) 

where 

a 0 = p^eW-p^eW, 
ai = 2k{p^E^U 2 -p (1) E^U 1 }, 
bi = a{E^-E^}, 
a 2 = k 2 {p^E^U 2 - p^E^U 2 } 
-^R 2 k 2 + m 2 -l), 

b 2 = ak{E^U 2 -E^U 1 }, 
where for the simplicity of notation, we used 
E$ = E$(k,R), (j = 1,2) 

where Em (k, R),(j = 1,2) are explained by 
(3.4)-(3.5). (i) When a = 0, (2.14) reduces to 


It is clear from the above inequality that the 
streaming has a destabilizing effect on the stabil- 
ity of a cylindrical interface, because Em is al¬ 
ways negative from the properties of Bessel func¬ 
tions. (ii) when a ^ 0 , we find that necessary 
and sufficient stability conditions for (2.14) are 

[3] 

61 > 0 , (2.18) 

and 

a^b\ — 016162 + 0261 < 0, (2-19) 

since ao is always positive . 0 

Putting the values of 00 , 01 , 02,61 and 62 
from(2.14) into(2.18) and( 2.19) we notice that 
the condition (2.18) is trivially satisfied since a 

is always positive, and from properties of Bessel 
( 2 ) 

functions Em is always negative. From (2.19), it 
can be shown that the condition for the stability 
of the system is 


R 2 k 2 +m 2 -l)+k 2 
R z 


pWpWE$E%\u 2 - Ur ) 2 

P&E$ - p^E^ 


n 

CLqLJ + CL\(jO + ( 2-2 = 0. 
Therefore the system is stable if 
(l\ — 4oq02 > 0, 


^(R 2 k 2 + m 2 - 1) 
R- 

2 p^p^E^E ( ^{U2-Ui) 2 
pV>E$ - p^E (2) 


(2.15) 

(2.16) 


E$E%\pV> -p^) 2 ' 
(E$ - E {2) ) 2 pWp( 2 ). 


( 2 . 20 ) 


The stability condition (2.20) differs from (2.17) 
by the additional last term: 


4MV 1 ) - p^ 2) ) 2 /[p W p^(E^ - Em) 2 ]- 
Thus the condition (2.20) is valid for infinites¬ 
imal a and when a = 0 the last term is absent. 


(2.17) 


We now employ multiscale expansion near the 
critical wave number. The critical wave number 
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is attained when a 2 = 62 = 0. The corresponding 
critical frequency, lo c is zero for this case. 

Introducing e as a small parameter, we as¬ 
sume the following expansion of the variables: 

3 

V = '}2( n r] n (9,z,t 0 ,t 1 ,t2) + 0(e 4 ), (2.21) 

n =1 
3 

4> {j) = '^2e n ^\r,e,z,to,t 1 ,t 2 )+0(e 4 ),(j = 1,2) 

n=0 

( 2 . 22 ) 

where t n = e n t(n = 0,1,2).° The quantities 
appearing in the field equations (2.3) and the 
boundary conditions (2.6), (2.11), and (2.13) can 
now be expressed in Maclaurin series expansion 
around r = R. Then, we use (2.21), and (2.22) 
and equate the coefficients of equal power series 
in e to obtain the linear and the successive non¬ 
linear partial differential equations of various or¬ 
ders. 

To solve these equations in the neighborhood 
of the linear critical wave number k c , because of 
the nonlinear effect, we assume that the critical 
wave number is shifted to 

k = k c + e 2 |U. 

3. First Order Solutions. 

We take 

$ ] = U jZ . (j = 1,2) 


The first order solutions will reproduce the lin¬ 
ear wave solutions for the critical case and the 
solutions of (2.3) subject to boundary conditions 
yield 

Vi = (3.1) 


AA - ( 

a 

01 ~ ( 

,P (1) 

AA - ( 

a 

01 ~ ( 


where 



r ) = 

E%\k, 

r) = 


+ ikU^j A(ti,t 2 )E i A ) (k, r)e 81? + c.c., 

(3.2) 

+ ikU-^j A(ti,t 2 )E^}(k, r)e J1? + c.c., 

(3.3) 


Im{kr)K' m (ka) - I' m {ka)K m {kr) 

I'm{kR)K' m {ka) - I' m {ka)K' m (kR) ’ 

(3.4) 

Im(kr)K > m (kb) - I' m (kb)K m (kr) 

I’ m {kR)K' m {kb) - I' m {kb)K' m (kR) ’ 

(3.5) 


d 

9 = kz + mO , l' m {ka) = — I m {kr) | r=Q , etc. 

with I m and K m are the modified Bessel func¬ 
tions of the first and second kinds, respectively. 


4. Second order solutions. 


With the use of the first order solutions , we 
obtained the equations for the second order prob¬ 
lem 


V 2 #=0, a = 1,2) (4.1) 

and the boundary conditions at r = R. 


[ d ^ ] dm 


\ dr dz 3 


- -w —Uj > -ar ]2 = 


,(j) 


4-2U 2 + ^+ 


R 


R 2 


0 . 0.2 


'A 2 e™ + P U) ^e^ + c.c. + 2a(^ + a 2 ) |4l| 2 , (j = 1,2) (4.2 - 4.3) 

01,7 m„ M 1 ' , , 

” U2 ^r- p Ul ^ + a \-a^ + 


' d 2 rj2 , 1 d 2 r]2 ( ?? 2 
R 2 AW 2 + R 2 


m 


p[— + ikU ) <J -1 - ( jp + ^ ) E m } + SaUki - 2 pU 2 k 2 
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a 


+ -^3 (R 2 k 2 + 2- 7m 2 ) \A z e 2lit + 


t ( — + ikU ) E n 


k\p 


8A 

dti 


_j_ cc 


a 


P[ —9 + k 2 U~ ) <! — 1 + I —y + k 


P* 


m 


R 2 


a 


+ m _ 2 )}>|A| 


The non secularity condition for the existence of the uniformly valid solution is 0 

dA 


dt\ 


= 0 . 


Equations (4.1) to (4.4) furnish the second order solutions: 

V 2 = -2 |A| 2 + A 2 A 2 e 2i ' 9 + A 2 A 2 e ~ 2i ^, 

</>? = B 2 ) A 2 e 2l ' } E^ l (2k,r) + c.c. + b^\to,ti,t 2 ), ( j = 1,2) 


where 


Ao = — 


— -<f 

11 (0, 2k, 2m) \ [ 


—pi2kU E 2m (5 + 2 (j& + + 1 } ( p + iW _) 


+^(2 + R 2 k 2 -7m 2 )}, 


+2 p(kU) 2 — i3akU 

B^ = ^ + l-^ + 2ikU j \A 2 , 

/3(J) = { jf) + iku > }{li- 2E ™ ) {je + fc2 ) } + 


«Q 2 


P 


U) 


(2 )db( 2 ' ) m <%W 


Otn 


wher eE%l = E^(2k,R). 


dtn 


/ a 2 

T \P 


p(^ + k 2 U 2 )\l-El(k,R)[^ + E 


— ~3 (k 2 R 2 + m 2 - 4 - 2ito 2 ) | |A| 2 , 


5. Third order solutions 

We examine now the third order problem: 

V^ } = 0 . (* = 1 , 2 ) 


(4.4) 

(4.5) 

(4.6) 

(4.7) 


(4.8) 

(4.9) 

(4.10) 

(4.11) 


(5.1) 


On substituting the values of rji, from (3.1)-(3.3) and rj 2 , 4> 2 from (4.6)-(4.7) into (A.7), we 
obtain 

4>\P = E^^JJk, r)A 2 Ae 1 ^ + E^\k,r)^-e 1 ^ + c.c., (5.2) 

where 


ri{j) _ 

°3 — 




x( — + a2 )+~‘(fc + 


1 f 2 2 + m 2 Em ( 3 m 


R 2 


R V R 2 


+ k 2 | [> [ ^ + ikU 


U) 
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f n 

\R 2 


m " 

~R 2 ~ k ^ 


+ ~YJ\ \ 4a 2 ( 4 + a 2 I - 3a 3 


P 


U) 


R 


-■ 7 "tt — ikUj 
o(3) J 


E^f^ + k 2 


1 

+ R 




(3 = 1 , 2 ) 


(5.3) 


We substitute the first- and second-order solutions into the third order equation. In order to 
avoid nonuniformity of the expansion, we again impose the condition that secular terms vanish. 
Then from (A. 8 ), we find 0 


where 


dD(0,k,m)dA 
i - + < 2crk c p + 


du 


dU 


pi r ^y~ + ikUjE m 


k c i/j, > A + qA 2 A = 0, 


+ 


+B 2 

3 o? 
R-p 2 


p(ikUC 3 E m + A 2 (^i~3kU - k 2 U 2 ^j 
(^ - iUk^j S.2E m E 2m (^ + k 2 ) - 1 } -*(^ + 2 « 2 ^ku(j + ikuj 


— E n 


m 

R 2 


+ F 


kU.(_a 


i ( 7—b 5 kUi ) + 


771 v 2 


-^{(2A 2 R-4:-4Ra2){l-m 2 )-2A 2 R(m 2 + k 2 R 2 )-^(m 2 + k 2 R 2 ) 2 +l(9m 2 + k 2 R 2 -6)}. (5.5) 


R 2 

3 

2 


Ei 


P 


+ 3k 2 U 2 - i2-kU 
P 


(5.4) 


We rewrite (5.4) as 
BA 

——b (ai + U2|A |~)A = 0, (5-6) 

ot 2 

which can be easily integrated as 

\A(t 2 )\ 2 = a lr \A 0 \ 2 exp(-2ai r t) 
x [air + « 2 r| Ao|' ® 2 r | A() |^ exp( 2fli r t)] , 

(5.7) 

where Ao is the initial amplitude and aj r = 
Maj, (j = 1,2) . 

With a finite initial value |Ao|, |A| may be¬ 
come infinite when the denominator in (5.7) van¬ 
ishes. Otherwise, |A| will be asymptotically 
bounded. The situation can be summarized as 
follows: 

(1) a 2 r > 0; stable. 

(i) a\ r > 0; \A\ 2 —> 0, as t 2 —► oo 

(ii) ai r < 0; \A\ 2 —> —ai r /a 2r , as t 2 —► oo 

(2) If a 2r < 0, 

(i) air < 0 ; unstable. 

(ii) air > 0, and |Aq| 2 > ~a\ r /a 2 r '■ unstable. 


(iii) air > 0 , and |Ao | 2 < —a\ r /a 2 r ■ stable 
and | A | 2 —> 0 as t 2 —» oo. Thus, a sufficient con¬ 
dition for stability is a 2r > 0 , which is due to the 
finite amplitude effect. The cylindrical system 
is nonlinearly stabile if ai r > 0 and the initial 
amplitude is sufficiently small. 


6 .Viscous asymmetric linear cylindrical 
flow 

In this section we consider the viscous poten¬ 
tial flow. For the viscous fluid, (2.12) is now 
replaced by 

pW(V^ (1) • VF) + V<^ (1) • \7F^J 

= / 2 ) (V^ (2) • VF) (^- + Vy? (2) • VF^j 

+{P 2 ~ Pi ~ 2/U 2 n • V <g> • n 

+2^in • V <g) V<£>^ • n + crV • n)|VF| 2 , (6.1) 
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where pi,p 2 are viscosities of fluid ’1’ and ’2’, 
respectively and we modify (2.13) accordingly. 
The nonlinear analysis for the viscous fluid is 
too onerous when the perturbation is asymmet¬ 
ric , we are content here with the linear analysis. 
Then linearizing (2.6), (2.11) and (6.1) we have 


0 





= 0 , 


(D ( d^ (1) dr/ dr] \ 

p v dr dt dzU ) 


ar], 


( 6 . 2 ) 


(6.3) 


P 


dt dz 


+ 2/x 


d 2 ^ 

dr 2 


(d 2 i] r] 1 d 2 rj\ 

l^cfe 2 R? r 2 d9 2 ) 


(6.4) 


When the interface is perturbed to r] = 
A exp [i(kz + mO — cot ), we recover the first order 
solutions (3.1)-(3.3), and the dispersion relation 


for the viscous fluid is same as (2.14), however 
a 0 = p^E^-p^E^, 


a 1 
bi 
02 


2 k{pWE$U 2 - p^E^Ud, 
a{E^ - E%>} + 2(/ii E™ - p 2 E^ 2) ), 
k 2 {p^E$U 2 - p^E^U 2 } 


a 


~^{R 2 k 2 + m 2 - 1 ) 
R- 


—2a 


pi E m _ 2 ) 


E 


62 = 


with 


* pi 2 ) 

-2fc(/x 1 C/ 1J E7 t (1) - M 2C/ 2J E7 t (2) ) 




and necessary and sufficient stability conditions 
are 

b 1 > 0, (6.5) 

and 

a 0 b'i - a 1 b 1 b 2 + a 2 bj < 0, (6.6) 

since ao is always positive. 


7. Numerical examples 

In this section we do numerical works using the expressions presented in previous sections for the 
film boiling conditions. The vapor and liquid are identified with phase 1 and phase 2, respectively, 



FIGURE 1. The critical wave number for m=l. 
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SO that T\ > Tq > T 2 . 

In the film boiling, the liquid-vapor interface is of saturation condition and the temperature To 
is set equal to the saturation temperature. The properties of both phases are determined from this 
condition. First, in figure 1 we display critical wave number k c , i.e., the value for which co = 0 in 
(2.14) Here we chose pi = 0.001gm/cm 3 , p 2 = 1 gm/cm 3 ,a = 72.3 dyne/cm, b = 2cm, a = lcm, R = 
1.2 cm, a = 0.1 gm/cm 3 s 0 



FIGURE 2. The stability diagram for the flow when m=l. The system is stable in the region 

between the two upper and lower curves. 



U2(U1=1) 


Fig.3.Viscous cylindrical flow for m=0.The region above the curve is stable region. 

From this figure we can notice that critical wave number increases as the velocity of fluid increases, 
the increment rate of the inviscid fluid being sharper at higher fluid velocities. In figure 2 we display 
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the region of stability of fluid in the nonlinear analysis as the velocity of one fluid increases while 
that of the other fluid remains unchanged. In these figures , u\ remains constant as 1 cm/sec while 
U 2 varies from 1 cm/sec to lOcm/sec. The region between the two curves is the region of stability, 
while in the region above the upper curve, the fluid is unstable. 0 

In Fig.3 and Fig.4 we present the results for viscous cylindrical linear flow.Here we chose 
pi = 0.0001gm/cm 3 , p 2 = lgm/cm 3 ,a = 72.3 dyne/cm, b = 2cm, a = lcm, R = 1.2 cm, a = 
.Igm/cm 3 s, p,\ = O.OOOOlpoise, p 2 = 0.01 poise 



Fig.4.Viscous cylindrical flow for m=l.The region above the curve is stable region. 


8. Conclusions. 

The stability of liquids in a cylindrical flow when there is mass and heat transfer across the 
interface which depicts the film boiling is studied. Using the method of multiple time scales, a first 
order nonlinear differential equation describing the evolution of nonlinear waves is obtained.With 
the linear theory the region of stability is the whole plane above a curve like in Fig.3,4, however 
with the nonlinear theory it is in the form of a band as shown in Fig.2. Unlike linear theory, 
with nonlinear theory, it is evident that the mass and heat transfer plays an important role in the 
stability of fluid, in a situation like film boiling. 

Appendix 

The interfacial conditions are given on r = R as 
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